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PERIODIC  MAGNETODIELECTRIC  LAYER  BOUNDED  BY  TWO 
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ABSTRACT 

The  problem  of  electromagnetic  wave  propagation  in  a doubly-periodic 
magnetodielectric  layer  bounded  by  two  uniform  infinite  media  is  solved  by  new 
method  based  on  the  rigorous  volume  integro-differential  equations  of  electromagnetics. 


The  Galerkin  method  is  applied  to  reduce  these  equations  to  a set  of  second-order 
differential  ones  with  constant  coefficients  in  field  functionals  which  contain 
information  about  geometry  of  the  scattering  structure.  The  special  scheme  of  equation 
set  solving  is  introduced  in  the  case  of  thick  layers  to  overcome  usual  numerical 
difficulties  associated  with  the  undesired  exponential  functions  in  the  expressions.  This 
method  unifies  the  treatment  of  both  TE-  and  TM-waves  by  replacing  £ by  p,  p by  e,  E- 
components  by  H-components,  H-components  by  E-components. 

METHOD 


Formulation  of  the  problem  is  as  follows:  from  the  region  2 ( z<-h ) with  complex 
relative  permittivity  £2  and  permeability  p2,  a linearly-polarized  plane  electromagnetic 
wave  is  incident  at  an  arbitrary  angle  (p  on  the  double-periodic  infinite  layer  (region  1) 
bordering  the  region  3 (z>0  ) with  complex  relative  permittivity  £3  and  permeability  p3 
(Fig-1). 


REGION  3 


P cv 


Figure  1.  Geometry  of  the  problem. 


Kii;v,  Ukraine,  IX-th  Internationa!.  Conference  on  Mathematical  Methods  in  Electromagnetic  Theory 


402 


MMET*  02  Pro(  ■/:/-:dings 


The  periodic  cell  of  the  layer  is  an  oblique-angle  parallelepiped  of  arbitrary  sizes  along 
the  r|i  axis  and  the  rj2  axis,  L,  and  Ln  are  periodicities  of  the  layer  in  the  rp  and  Tq2 

direction,  respectively.  The  parallelepiped  is  characterized  by  the  complex  relative 
permittivity  £|  (r| i , p2)and  permeability  pi  (pi,  r\2)  and  has  the  thickness  h.  We  suppose 
that  the  incident  wave  is  TE-polarized  and  a is  the  angle  between  the  r\\  axis  and  the 
electric  field  vector  £ lying  in  the  plane  of  the  layer  and  consider  the  field  components 
in  the  (p,  z)  ortogonal  coordinate  system  connected  with  the  incident  field 
polarization.  We  have  to  obtain  the  transmitted  and  reflected  fields  in  the  immediate 
vicinity  of  the  layer.  The  form  of  the  integro-differential  equations  for  the 
electromagnetic  field  is  taken  from  [1]: 


£(r)  = £0(r)  + 1 (VV  + k~)  f(^'  -1  )E(r')G(r-r')dr'  + 
4rc  £, 

+ lk]  h-Vx  f(  h -\)H (r')G(r-r')dr\ 

4n  V £ , p, 

fHr)  = n0(r)+  - (VV  + ^)f(^  -\)Ft(r')G(r-r')clr'~ 

4 n p, 

_ Lki  . v x r')o(\r  - r' )dr' . 

4n  p,  £, 


Here  the  Green’s  function  is  presented  in  the  integral  form,  V is  the  scatterer  volume. 
kt  =k-sj e2p2  and  £ is  a wave  number.  For  beginning  we  suppose  that  the  periods  are 
partitioned  into  segments  with  constant  material  parameters.  For  each  period  segment 
numbered  (kj)  we  can  write  the  notations  b^'  <p'A  < 6* , b'}  ' <pi  < b^  , £,(>//)  = £w 
and  p,  (>/, ) = pw  • 

Following  the  algorithm  given  in  [2]  we  present  the  field  in  each  period  segment  as  an 
expansion  in  terms  of  the  spatial  harmonics  numbered  (r,s)  and  act  on  the  equation  for 
these  fields  by  the  linear  operator 


Ak!  F(y\'vr\'-,,z')  - 
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Thus  we  can  obtain  the  set  of  linear  differential  equations  for  the  field  functionals.  By 
summing  these  equations  over  all  segments  we  can  express  the  fields  and  field 
functionals  for  individual  segment  through  the  ones  for  another  segment,  e.g.: 


rW  _ CrCti  1) 


( A A, , )(^p->x,y, 


r A 1 £■’. 


kT 
1 V)(l  ' 


where 
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Next  we  solve  the  equation  set  for  the  field  functionals  and  thus  for  the  field 
components  (the  procedure  in  paper  [2]).  Using  the  extinction  theorem  at  the  last  stage 
we  express  the  internal  fields  of  the  structure  (region  1 and  3)  through  the  incident  wave 
field.  Upon  separation  the  equation  set  to  be  solved  onto  the  subsets  including  the 
exponential  functions  with  a positive  and  negative  real  part  of  eigenvalues,  we  apply  the 
iterative  procedure  developed  for  calculations  of  characteristics  of  scattering  from  a 
very  thick  losser  doubly-periodic  magnetodielectric  layer  in  a wide  frequency  range. 

RESULTS 

Presented  in  Fig. 2 are  the  results  of  the  transmission  coefficient  calculation  for  the 
periodic  structure  composed  of  the  array  of  square  parallelepipeds  with  e = 3+0.0  li, 
p = 1,  lying  on  the  half-space  with  £3  = 3,  p.3  = \,h/ =10,  = 0.5 

IT)i  = Z,T]i  = Lt] , a = 25°, (p  = 0.001°.  The  example  demonstrates  the  iterative  procedure 
stability  for  very  thick  doubly-periodic  structures. 


Figure  2.  Transmission  of  the  normally  incident  plane  wave  through  the  thick  doubly- 
periodic  array  lying  on  the  half-space. 

Our  results  are  in  a good  agreement  with  data  concerned  thick  periodic  structures  and 
presented  by  other  authors  (see,  for  example,  paper  [3],  where  the  one-periodic  grating 
of  rectangular  rods  with  h! L~  4 was  considered). 
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